load Z_a.dat;
load PHI.dat;
load RHO.dat;
load U.dat;
load V.dat;
load P.dat;
%U=U(end:-1:1,:);
%V=V(end:-1:1,:);
line = 600;
xa=linspace(0,15,line);
ya=xa;
[x,y]=meshgrid(xa,ya);
figure(1);
mesh(x,y,Z_a);
figure(2);
mesh(x,y,PHI);
figure(3);
mesh(x,y,RHO);
figure(4);
mesh(x,y,P);
figure(5);
quiver(x,y,U,V);
